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^ ; Abstract 

We investigate quantitatively the extent to which having a primordial 
(Sj ! leptonic asymmetry (n v ^ rip) relaxes the bounds on light sterile neu- 

trinos imposed by BBN and LSS. We adopt a few assumptions that 
\q ; allow us to solve the neutrino evolution equations over a broad range 

of mixing parameters and asymmetries. For the general cases of ster- 
43 ■ ile mixing with the electron or muon neutrino, we identify the regions 

that can be reopened. For the particular case of a LSND-like sterile 
neutrino, soon to be rejected or confirmed by MiniBooNE, we find that 
an asymmetry of the order of 10~ 4 is needed to lift the conflicts with 
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cosmology. 



^ : 1 Introduction 



Cosmology currently poses some of the most stringent bounds on the properties of sterile neu- 
trinos v s - Standard Model (SM) gauge group singlet fermions - namely, their mass and mixing 
parameters with their active counterparts v e , and v T . 

In particular, by comparing the predictions of the yield of primordial light elements with 
observations, Big Bang Nucleosynthesis (BBN) efficiently constrains the possible manifestations 
of sterile neutrinos. Having a greater than 3 effective number of neutrinos would increase the 
rate of expansion of the universe and leads to larger 4 He and Deuterium abundances. Even if 
steriles are produced only after the decoupling of the active neutrinos from the cosmological 
plasma, so that the effective number of neutrinos remains roughly constant at 3, they would 
lead to a depletion of electron neutrinos and anti-neutrinos, which also boosts the amount of 
4 He and Deuterium produced P . In addition, observations of the matter power spectrum of the 
Large Scale Structures (LSS) in the Universe limits the total mass in neutrinos today, including 
the sterile ones 0. 



However, naturally light or cleverly lightened sterile fermions arise in many theories beyond 
the Standard Model. On the concrete experimental side, moreover, the result of the LSND 
experiment [3] still stands. The most common interpretation requires the existence of a sterile 
neutrino with a mass (9(1 eV) mixing with both the electron and muon neutrino with angles 
O(10~ 3 ). Such parameters are solidly ruled out by the constraints of standard BBN and by 

lss HEHaE]. 

As well known, the same region of the parameters space is currently being scrutinized by 
the MiniBooNE experiment jH], which is expected to deliver results soon. Should MiniBooNE 
confirm the existence of the LSND-like sterile neutrino, this would surely signal the need of 
a drastic revision of standard cosmology. If instead MiniBooNE strengthens the experimental 
exclusion contour beyond the LSND region, the bound from BBN on this particular mixing 
configuration is still stronger - a few orders of magnitude in mass and mixing angle - and these 
issues remain: in what ways and how much of the parameter space can be reopened? 1 

The bounds from cosmology are based on the assumption of "standard evolution". They 
relax considerably if some non standard modification is postulated. One such modification that 
has attracted much attention is to invoke the presence of a large primordial asymmetry for one 
or more of the active neutrinos seas. Of course, one might expect the value for a primordial 
asymmetry to be set by that of the baryonic asymmetry 77 = (tib — n^/n^ ~ 6 • 10~ 10 . To respect 
charge neutrality, the asymmetry in the charged leptons must be of the same order. But there is 
no such constraint on the asymmetries of the electrically neutral neutrinos, and they are therefore 
allowed to take values that can be orders of magnitude larger, up to 0(1CT 2 ) [TU]. Indeed, there 
are many examples of cosmological scenarios that lead to a large primordial neutrino asymmetry 
while preserving the observed smallness of the baryonic one, such as GUT scenarios, the Affleck- 
Dine mechanism, Q-balls, etc. [H]. From the purely phenomenological point of view, this is thus 
an important possibility and it has relevant consequences on the sterile neutrino bounds. 

The basic idea is the following [T2j. A non vanishing asymmetry L v produces an additional 
term in the matter potential for the active neutrinos that is proportional to the asymmetry 
L v itself. This increases the magnitude of the diagonal entries of the real part of the effective 
Hamiltonian written in the flavor basis relative to its off diagonal ones, thus suppressing the 
active neutrino oscillations into the sterile states. Since sterile neutrinos can only be created 
through active-sterile oscillations, this means they will be less efficiently or not produced at all 
in the Early Universe. 

The main goal of this paper is to study quantitatively the extent to which such a mechanism 
works. Some aspects of this issue have been addressed before (see e.g. O EH EH EH])- We 
aim at an analysis that includes as much content of the full neutrino evolution equations as 
possible, such as the established active-active neutrino mixing, and that covers a broad range of 
parameters. We compute the effects in terms of observable cosmological quantities, such as the 
yield of the primordial elements, and compare them with the most recent results, which have 
gone through significant changes. 

The paper is organized as follows. In SecEl we describe the neutrino mixing formalism, 

1 This could also be of importance for the r— process nucleosynthesis in SuperNovae, as Ref. [3j finds that the 
region of successful synthesis in the presence of an active/sterile oscillation has some overlap with the LSND 
area. 
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identify the parameter spaces we are interested in, and present the evolution equations that we 
solve. In SecOJwe present our results and in Sec|3]we summarize our conclusions. Appendix 1X1 
lays down some relatively standard details that govern the evolution of the temperature of the 
primordial plasma in the time of interest. 

2 Formalism and computations 

In this Section we delineate the neutrino mixing formalism and the computational strategy. The 
uninterested reader can skip to SecOJ the basic points are that, for any given set of values 
of the active-sterile mixing parameters (Amf 4 , 9 S ) and primordial leptonic asymmetry L° = 
(n° — n°)/n° (taken to be the same for the z/ e , and v T seas, for the sake of defmiteness 2 ), we 
compute the evolution of the densities of the four species of neutrinos and antineutrinos, from 
an initial high temperature down to a fraction of an MeV. With these, we follow the evolution 
of the neutron to proton ratio n/p: it determines the yield of 4 He and Deuterium, which we 
compare with observations. On top of that, the total amount of produced sterile neutrinos is 
subjected to the bounds from LSS. 

2.1 Neutrino mixing 

To parameterize the mixing of the sterile neutrino with one specific linear combination of the 

3 active states, we adopt the general formalism introduced in Ref. [I]. The extra (mainly 
sterile) neutrino has a mass m 4 . It is allowed to mix, with a mixing angle 8 S , with an arbitrary 
combination of the flavor eigenstates of the active neutrinos, represented by a complex unit 
3- vector n 



The 4x4 neutrino mixing matrix V, that relates flavor to mass eigenstates as ^ e ,/i,r,s = V • ^1,2,3,4 
is expressed in this parameterization by 



where U denotes the usual 3x3 mixing matrix: Vf t = Uuu mv with t = e,jj,,T and / and 
m denoting flavor and mass eigenstates respectively. It contains the mixing angles #12, #23? $13 
among the active species. Their values, together with those of the active-active squared mass 
differences Am^ 2 , jAm^l, are fixed by the solar, atmospheric and terrestrial neutrino experiments 
(see [Tj\ for a recent compilation). For simplicity, we restrict to the case of normal hierarchy 
(A777J3 > 0) and assume #13 = 0. Since oscillations are sensitive only to mass differences, without 
loss of generality we set m\ = 0. This is also the most conservative choice with respect to the 
bounds from LSS. The oscillation results will be presented in terms of the active-sterile squared 
mass difference Am^. 

2 In principle, different primordial asymmetries for the different flavors and even conspiracies that compensate 
the net effect arc possible. Our choice seems more generic and it is made on the basis of a balanced initial 
condition among flavors at high energies. The equalization discussed in |1(J| will also eventually occur, albeit it 
could take place at low temperatures. 



n ■ V = n e u e + n^v^ + n T u T . 



(1) 




(2) 
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In the following we will focus on the particular cases of sterile neutrino mixing with a pure 
electron neutrino state ("es mixing") and with a pure muon neutrino state ("/is mixing"). These 
correspond to ng = (1, 0, 0) and n t = (0, 1, 0) respectively. We will not address the case of mixing 
with a pure tau neutrino because it is equivalent to the muon case, as they experience the same 
matter potential and have the same roles in weak interactions. We will limit our investigation 
to Am| 4 > 10~ 4 eV 2 ~ Am 2 2 for the es mixing case and to Am 2 4 > 10~ 2 eV 2 ~ Am^ for the 
/is mixing case so that the (mainly) sterile state is always heavier than the flavor with which 
it mixes (often called, in absence of asymmetries, the "non- resonant case"). Finally, we will 
consider the case of the LSND sterile neutrino. This case requires a mixing with both v e and 
z/ M , and therefore interpolates between the previous two cases we have discussed. 

2.2 Neutrino evolution and BBN computations 

The machinery of Big Bang Nucleosynthesis predicts the abundances of the primordial elements 
produced by the synthesis that occurs at T B bn ~ 0{1 to 0.1) MeV. We will focus on 4 Helium and 
Deuterium. The resulting abundances are dependent on the baryon asymmetry 7] = {n-Q— n B )/n 7 
and on the density of relativistic species in the early Universe - photons, and neutrinos. Fixing 
the value of rj with CMB data 3 to be rj = (6.11 ± 0.27) ■ 10 -10 (we adopt the value from pj, fully 
compatible with that from [IB]) and taking the energy spectrum of photons and to be those of 
the standard Bose-Einstein and Fermi-Dirac ones throughout the entire period of evolution, we 
need to compute the evolution of the densities of the four species of neutrinos and antineutrinos 
during the whole period of BBN, from T in T BBN to T fin < T B bn- This evolution is shaped 
by several processes occurring to the neutrinos: the active species oscillate among themselves as 
dictated by the active-active mixing angles and mass differences and by the matter potentials of 
the primordial plasma; the sterile species are populated by oscillations from the active species 
as dictated by 8 S and Am 2 4 ; eventually active neutrinos decouple from the thermal bath, just 
before e + e~ annihilate into photons. In the rest of this Section we delineate the computational 
formalism that allows us to follow this evolution. 

2.2.1 Treatment of the neutrino evolution equations 

The most convenient tool to follow the system of neutrinos is that of the density matrix [TH] . 
Intuitively, its diagonal entries correspond to the populations of different species of neutrinos, 
while the off-diagonal entries express the quantum superposition of flavors. Written in the 
flavor basis, the evolution equations for the full 4x4 neutrino matrix g(p, t) and antineutrino 
matrix g(p, t) with momentum p take into account (i) the vacuum oscillations (active-active and 
active-sterile), (ii) the matter effects in the primordial plasma, including those due to the lepton 
asymmetry that we are introducing, (iii) the ve <-» ve scattering reactions and the vv <-> e + e~ 

3 We are assuming that the neutrino asymmetries that we are introducing have a too small effect on CMB to 
be able to modify this value, as can be inferred from the analysis in |23) . 
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annihilation reactions. They read 1^ 122] 



scattering & 
annihilation, 

scattering & 
annihilation, j2J 



The Hubble parameter 



H 



8n 3 G 



N 



"Ptot 



(3a) 
(3b) 

(4) 



depends on all the energy densities p, including those of active and sterile neutrinos and 
antineutrinos, which are functions of the diagonal entries of the density matrices: p to t = 
p 7 + p e ± + Yn=i PniPii) + Yn=i PniPa)- Tne Hamiltonian in matter 



H = Ho + diag(V e ,V M ,V r ,0) 
is composed by the vacuum Hamiltonian in the flavor basis 

H = —VM 2 V ] M 2 = diag (ml ml ml ml) 
2p 

and the matter potentials V; for each flavor 

"1 



(5) 



(6) 



V e = ±V2G F n 1 

8V2G F 
3M 2 , 



-77 + 2L Ve + L Vfi + L Vr 



p({E e )n e + (E s )n s ) 



8\/2Gp / , . , . N 

3M 2 V[\Eu)n Ve + {E p )n Pe ) 
z 



V M = ±y/2G F n 1 
V T = ±V2G F n 1 

v., = o 



V + U e + 2L V/l + L Vt 

V + L Ue + L v + 2L Ut 



8 \J~2G p i . . / 7-i \ \ 

3M 2 P{\ E ^ n ^ + (^)^J 

8^/2G p / . . . . x 

3A//2 P{{Ev)n UT + (£p)n p J 



(7a) 

(7b) 

(7c) 
(7d) 



where we are making use of the standard definitions of the asymmetry 



L v 



and of the number density and average energy of a species x with distribution f x (p) and g x spin 
degrees of freedom 



n 3 



/d?p 1 f d?p 

frtf 9 * Uvl {Ex) = nJWf 



g x E(p) f x (p). 



(9) 



For photons, E(p) = p and n 7 = 2((3)/tt 2 T 3 . Neutrinos and antineutrinos are distinct and 
have g v — g v — 1; also, because they are nearly massless relative to the ambient temperature, 
E{p) = a/ ml + p 2 ~ p. 



The first terms in each expression of the V; are the contribution from the asymmetries: the 
baryon asymmetry rj and the asymmetries in the background neutrinos of all flavors, with a factor 
2 in front of the same-flavor contribution; the plus sign applies to test neutrinos, the minus sign 
to test antineutrinos. Notice that the leptonic asymmetry will start from its primordial value L°, 
one of the free parameters of our analysis, but it is in principle a dynamical quantity, determined 
by the evolution of the neutrino and antineutrino densities. This might complicate the equations 
non trivially. We will further comment on this point below. 

The second term is the extra contribution due to the presence in the plasma of background 
neutrinos of the same flavor as the test neutrino. This term does not vanish even in case of 
symmetry in the background densities and effectively provides a thermal mass to neutrinos 
immersed in the bath. The third term in the electron neutrino potential similarly comes from 
the background electrons and positrons of the bath, and it is not there for neutrinos of other 
flavors. 

The last terms in eq.s Q contain the rates T's of the ve <-» ve, vv «-> vv scattering reactions 
and the vv «-> e + e~,vv annihilation reactions. Their complete form can be found in (25- We 
will provide the explicit form that they take in our formalism below. 

The momentum dependent equations eq. (J3J) do not allow for an efficient way of studying 
the evolution of neutrinos and antineutrinos over a broad range of parameters, which is the 
goal of our paper. We work therefore in the limit in which the neutrino distribution can be 
approximated as Fermi-Dirac functions with a negligibly small chemical potential, multiplied 
by an overall normalization factor that can be interpreted as an average of the available phase 
space: 

eta*)- 5^1- ^-i^TT ( 10 > 

Here p and p are 4x4 momentum-independent density matrices that evolve with time only. 
Even though we are dealing with particle asymmetries, neglecting the chemical potential is a 
good approximation as long as the asymmetries under consideration are small, as will be in this 
paper. Quantitatively, one can estimate a priori the typical size of the asymmetries that we will 
be considering by comparing the contributions of the asymmetry term and the thermal term in 
eq. 1)17]) . For a neutrino of average momentum p ~ 3.15 T, the two become comparable for 

L v « 3.15 ^ 77r L (l + -i-) « 3 1(T 7 ( — (11) 



180 C(3) V 1 + c 2 vJ Mf ~ 3 10 UoMevJ 



In other words, an asymmetry L u of the order of (9(10~ 6 ) is already enough to give a dominant 
contribution to the matter potential (and therefore an important effect on active-sterile oscil- 
lations) in the range of temperatures T ~ few tens MeV and range of masses Am^ < few eV 2 
in which we are interested, as one can verify ^3] that the production of sterile neutrinos (in 
absence of asymmetries) dominantly occurs at a temperature 

/nrp \ 1/3 /Arr) 2 \ 1 ^ 6 

T prod „. « lOMeV ly^j (cos2^) 1/6 ( (12) 



Since an asymmetry L v translates into a dimensionless chemical potential £ = p/T via 

12C(3) V T J V 7r 2 



= (13) 
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for small values such as 10~ 6 , we have £ ~ L u , which requires a negligible modification to the 
pure FD distribution. For much larger asymmetries our assumption would be questionable. 

Other kinds of deviations from a pure FD distribution could occur during the evolution, 
possibly due to the conversions into sterile neutrinos themselves. These have been shown e.g. 
to be not negligible for Am 2 < 10~ 8 eV 2 j2E], which is however a region far from that of our 
investigation. More recently, Ref.s ^H] have shown that matter resonances post weak decoupling 
might leave both active and sterile neutrinos with a highly non-thermal spectrum for some choices 
of neutrino parameters and energies. This appears to be the case in particular for values of the 
initial asymmetries larger than those that we will investigate. However, a fully momentum- 
dependent evolution analysis, which is beyond the scope of this paper, would be necessary to 
assess the importance of this effect in the entire range of asymmetries and mixing parameters 
that we consider. 

With the assumptions outlined above, the evolution equations for the neutrino and antineu- 
trino q matrices can be integrated over the momentum p and reduced to the corresponding 
equations for the p matrices. Explicitly, they read 

P=^% = ^\H iP \- \t« ( {Gl p-D- 2G S ( P - nc s 

+ {Gl (p - p^)} + 2G a (p - jf*)G^ , (14a) 

+ {Gl, (p - n } + 2G a (p - p cq )G a ) , (14b) 

where now 

« = «o + diag(V e ,V„,V x ,0), (15) 

w » = lwrm VM2v *' (16) 
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V e = ±V2G 



1 3 / T v \ ( 

" V + 8 v ~T ) \ 2<ypee ~ Pee ^ + ^ Pwt ~ p ^ + ^ p 



TT \ J TT 



196 tt 2 C(4) 
180 C(3) 

196 vr 2 C(4) 



V2d 



'Ml 



T 4 + -T^ COS 2 6» w (p ee + p ee ) 



4 



(17a) 



1 3/T, 

2 8 ( I 7 



(Pee - Pee) + 2(p w - p w ) + (p rr - p rT ) J 



V2G, Tv ^ 4 ~ 



(17b) 



180 C ( 3 ) v ^^m 2 / 4 T ' COs2 ^ (Pw + ^ ; 



V T = ±V2G 



1 3 fT u \ 3 ( 

2 V + 8 i Y ) V Pee ~ Pe ^ + ^ ~ P ^ + 2iypTT ~~ Pl 



196 tt 2 C(4) r- T v 1 4 2 
^80~C(3) F 772-7 " COS ^ PrT PrT ^' 



(17c) 



M^4^ 

sin^ = 0.2312 is the weak mixing angle. 
In these equations we have denoted 

r° = G F Tl (18) 

while G s and G a are diagonal matrices that contain the numerical coefficients for the scattering 
and annihilation processes for the different flavors. Explicitly 



G s>a = diag(7 s e a , nft ^ 0) with 



(7 s e ) 2 = 



3.06, 
= 2.22, 



(7a 6 ) 2 = 



0.50, 
= 0.28. 



(19) 



The numerical coefficients 7 can be determined from the expressions for the reaction rates 
reported e.g. in [2"4] . 

Our choice for the scattering and annihilation terms deserves some comment. They are 
based on eq.(j2J) integrated over all momenta. In particular, its form can be obtained by consid- 
ering the density matrix as a function of its difference from its equilibrium value, p = p eq + 5p 



and p 



cq 



6p, where p eq = diag(l + 4/3L„ eJ 1 + 4/3^, 1 + 4/3L„ T , 1 + 4/3L Vs ), p eq 



diag(l — 4/3L„ e , 1 — 4/3L„ M , 1 — 4/3L„ T , 1 — 4/31^). Strictly speaking, the equilibrium density 
matrices are the solutions that would make the RHS of eq.([3*| vanish, which means p eq and p eq 
must be proportional to the identity matrix: L Ue = L v = L Vt = L Us = L v . However, at high 
temperatures when oscillations are heavily suppressed 7i becomes nearly diagonal, and equilib- 
rium can be achieved with simply a diagonal matrix with possibly different asymmetries L v .. 
Furthermore, because the only occurences of p cq and p eq are those inside the anti-commutator 
of the annihilation terms 4 , and because (G Bia ) M = the elements (p cq ) ss and (p eq ) ss do not 
appear anywhere in the equations and hence can be set to zero. This is consistent with the 



^Because the matrix G s is diagonal, the scattering terms do not depend on p° q and p eq , as can verified by 



writing them as nested commutators (l/2)r° (— {G 2 S , p — p cq } + 2G s (p 
(l/2)T°[G s> [p,G s }}. 



p eq )G s ) = {l/2)T°{G s ,[p-p^,G s 
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fact that sterile neutrinos can be driven to equilibrium with other neutrino flavors only through 
oscillations. 

In general, the effect of the scattering and annihilation terms is to drive the density matrices to 
their equilibrium form, effectively replenishing the active neutrino populations at the expenses of 
the e + e~ bath and damping the flavor superpositions. When the rates T Uj . (that strongly depend 
on T) are overwhelmed by the expansion of the Universe (r„ < H) the so called neutrino freeze- 
out occurs, though neutrinos still take part in the n <-> p reactions (see below). It is not difficult 
to see that the diagonal entries in the scattering terms vanish, hence the fact that scatterings 
do not contribute to the change in the number density of neutrinos and antineutrinos is built in 
the equations. On the contrary, the off-diagonal entries are affected by the total scattering rate 
because both elastic and inelastic scattering play a role in damping out the coherence between 
different neutrino flavors. 

Moreover, we note that the equations by themselves guarantee the conservation of the total 
lepton number, given by the algebraic sum of the lepton asymmetries of the 3 active and 1 sterile 
flavors (which can individually vary). To see this, subtract the v from the v equations to get 

• - w i , Toy -i 1 rof [0° ~ P cq )G s ~ h.c.,G s ] + V terms + \ , . 

P - P = -i[H,p] + An, P ] - -r ^ [{p _ pCq)Ga _ h ^ GJ + _ termg j . (20) 

The commutator structure ensures that the total lepton number is strictly conserved: 

d 4 

Jt £>« - p«) = Trace ([- ■■>■■■]) = °- ( 21 ) 

i=l 

We finally stress that the popular (and simpler) anticommutator approximation of the scat- 
tering and annihilation terms (see e.g. eq.(283) of [22]) cannot be applied to our investigation of 
lepton asymmetries without additional constraint equations because it does not conserve total 
lepton number. 5 

One last comment is in order. Although the evolution equations presented in eq. (114)1 treat the 
asymmetries p aa — p aa as fully dynamical quantities, for numerical reasons we find it necessary 
to restrict the analysis to the case where the asymmetry is constant 6 . This is a reasonable 
approximation as long as the evolution itself does not induce drastic modifications. Indeed the 
latter have been studied in some detail and found to be important in the region of small Am;? 
and mixing angle [To^ I27j. but not in the main region of our interest. With our simplification, 
the evolution of neutrinos and antineutrinos are however not fully independent. In the event 
that a positive signal for a sterile neutrino is found at the MiniBooNE experiment or elsewhere, 
it might be worthwhile to abandon the above simplification and tackle the numerical problem 
of solving the evolution equations in their full complexity for the parameters possibly selected 
by that discovery. 

The determination of ^ is in principle quite involved and we leave its discussion for Ap- 
pendix [X] For the purpose of the above equations, we make use of the standard expression 

5 Whcn L v = 0, the difference between the anti-commutator approximation and our equations here can most 
iikely be attributed to our inclusion of Pauli blocking factors. 

6 As a technical aside: we implemented this by replacing p with p— (8/3)L in the commutator —i[H, p], where 
L = diag(L t/ , L v , L v , 0). By subtracting the p equations and these modified p ones, it is straightforward to show 
that the diagonal terms are zero. 
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T = —H{p Vtot ) T, where H contains the temperature dependent total energy density, including 
that of all neutrinos. The standard definition of p tot is recalled in eq. ()36a|) . In the density 
matrix formalism 

N v = 2^ a +p aa ) (22) 



We therefore solve these evolution equations from an initial temperature of T = 200 MeV, 
much higher than the typical production temperature of sterile neutrinos for any Am 2 in which 
we are interested (see eq. (JUJ) to a final temperature T <C 1 MeV. At 200 MeV, the abundance 
of the sterile neutrinos is assumed to be zero, for even if they are produced by some New Physics 
mechanism at very high scales, they would have been diluted by the reheating of the Universe 
due to the annihilation of the SM degrees of freedom). 7 The abundance of active neutrinos and 
antineutrinos at the initial temperature are 1 + (4/3) L u and 1 — (4/3) L u respectively. 



2.2.2 Neutron/proton evolution equations and primordial elements 

After solving the neutrino densities evolution as a function of temperature, we can study the 
relative neutron/proton abundance, which is the all-important quantity for primordial nucle- 
osynthesis, since they are the building blocks of the nuclei and essentially all neutrons are incor- 
porated into some light element in the process. The neutron abundance at the moment that the 
synthesis begins practically fixes the proportions of all the products, n/p evolves according to 

dT dv TL 
^^dT =Tp - {1 - r) - rT — r = ^f^ p (23) 

r p _+ n is the total rate for all the p — > n reactions and T n ^ p for the inverse processes. They 
depend on the electron neutrino and antineutrino densities computed above. Explicitly they are 
given by (see e.g. 



(1 + 3g 2 A )G 2 F /(2n 3 ) ( [ dE e p 2 u E e p e f De (p u ) [1 - f e+ (E e )]\ Eu=Ee+Am 

f°° (24) 

+ / &Pv PlPefe(E e ) [1 - f„ e (p„)]\E 

Jo 

/•Am 

+ / dE e plE e p e f Pe (p u )f e (E e )\ Pu=Am - Ee 

J m e 



7 For the opposite possibility, see |25| . 
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(/■Am 
/ dE e plE ePe [1 - UM\ [1 - f e (E e )}\ Pv=Am . Ee 
J m e 

r 2 ( 25 ) 

+ / d Pu P iE ePe f Ue ( Pu ) [1 - f e (E e )]\ E 

Jo 

poo 

+ / dE e plE ePe f e+ (E e ) [1 - UAPu)]\ p 



that we evaluate taking standard FD distributions with vanishing chemical potential for electrons 
and positrons, and the expressions in eq. fllOl) for the neutrino distributions. Here Am = m n —m p , 
the mass difference between neutron and proton, qa = —1.267 is the axial coupling for the 
neutron and the overall coefficient is best determined in terms of the standard neutron lifetime 
in absence of a thermal bath 



T ' n 2vr 3 



m 5 e \ , \ = J deeV^l\e- — \ ~ 1.633, (26) 



experimentally measured to be (885.7 ± 0.8) s [2 



The determination of ^ is affected by the several phenomena that go on in the range T ~ 
1 MeV which is under examination. In particular, on top of the expansion that cools the Universe, 
there is the production of possible extra degrees of freedom (the sterile neutrinos) and the heating 
of the thermal bath due to e + e~ annihilations. We deal with this in Appendix El 

From the final value of n/p a set of nuclear physics processes determine the predictions for 
the synthesis of the light elements. To perform this step we use the procedure briefly described 
in [1], calibrated to agree with more refined codes. The predictions are finally tested against 
the observations to draw constraints on the allowed modifications. We focus on the 4 He and 
on the Deuterium abundance: the observational determinations of both quantities are plagued 
by controversial systematic uncertainties that, especially in the case of 4 He, seem to restrict the 
margin for improvement; the case of Deuterium may have brighter prospects in this regard. We 
adopt the latest PDG recommended values jSO] for the observational determinations of: 



Y p = 0.249 ± 0.009, 
Y D /Y H = (2.78 ± 0.29) • 10" 



(27) 



where Yx = nx/n E and Y p is the traditional notation for Kt He . These values are dominated by 
the results in jHI] and [22] and are considered quite conservative (see [3H] for a slightly expanded 
account). For this reason we adopt 



Y p < 0.258, 



Y D /Y U < 3.07 ■ 10~ 5 
as current robust observational upper bounds. 



{21 
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2.3 Large Scale Structure bounds 

The measurements of the matter power spectrum of Large Scale Structures from a variety of 
techniques, combined with CMB data, allow one to constrain the amount of energy density in 
neutrinos today, including the sterile one (see e.g. [2]). The bound translates into a function 
of the mass and abundance of the sterile neutrinos [33]. In our case these two quantities are 
related, because the neutrino mass differences (together with the other parameters) determine 
the final abundance. If the sterile neutrino is fully populated by oscillations, its mass contributes 
fully to the computation of the total energy density. If it is not efficiently produced, the bound 
on its mass can be somewhat lifted. This is what we will study quantitatively below. 

In the regimes in which we are interested, the bound can be expressed to a good approxima- 
tion in terms of the quantity 

= 1/2 tt[Af ■(, + ;)] 

93.5 eV V ; 

where M is the 4-neutrino mass matrix introduced before and the neutrino densities p, p are 
taken at late cosmological time, in our case the end of the evolution. We assume here that its 
value is not significantly modified at low T. Here as usual h = Hq/{100 km/s Mpc). We adopt 
the constraint 

Q u h 2 < 0.8 ■ 10" 2 , (30) 
based on j3], which is in agreement with most of the recent analysis 



3 Results 

As discussed above, we solve numerically the neutrino and neutron/proton evolution equations in 
a broad range of values of the active-sterile mixing parameters and of the primordial asymmetry 
L v . We focus on the case of a negative value for the primordial asymmetry, motivated by the 
following. Because steriles are produced solely through oscillations from active states, their rate 
of production up to dimensionless factors goes as p ss ~ sin 2 (2# mat ter)r , where for small angles 
sin (2$ ma tter) ~ tan 

(2# ma tter) = sin 2 20 fl /(cos20 s - 2(£„/Am 2 )(V thcrma i + V asymm )) 2 [21 EE]. 
From eq. (JSJ) and eq. (jl7[) . a negative asymmetry would make all terms in the denominator of 
tan 2 (26 , mat tor) positive, leading to a smaller effective mixing angle and hence a reduced production 
rate. For large values of L u , the asymmetry term, which changes relative sign for v, could become 
dominant over the thermal term and allow for resonant transitions in the anti-neutrino channel 
to take place. As part of our approximations discussed above, we are neglecting these effects. 

We see noticeable suppression of v s production when the absolute value of the asymmetry is 
above C(10~ 7 ), consistent with the estimates in eq.flTTJ). The suppression causes the population 
of sterile neutrinos to rise later in the evolution, or to be kept at zero in the limiting case. 

3.1 Bounds on sterile mixing parameters 

Fig.sQshow the bounds on the planes (Am 2 4 , tan 2 6 S ) for the mixing cases es and ps, defined in 
Sec|2] We report the contours for Y p = 0.258 and Q u h 2 = 0.8 10~ 2 , so that the portions above 
and to the right of these contours are to be considered excluded according to eq. ([2~%]) and eq. (|3*U|) . 
We find that the bound from Deuterium abundance is less constraining than that from 4 He, so 
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Figure 1: The bounds from BBN and LSS on sterile neutrinos mixing with electron neutrinos (left 
panel) and with muon neutrinos (right panel), for different values of the primordial asymmetry L v . 
The regions above and to the right of the thin lines are excluded by BBN because they correspond to 
Y p > 0.258. The regions above and to the right of the thick lines are excluded by LSS because they 
correspond to Q v > 0.8 10 -2 . The gray dotted lines enclose the excluded regions from a combination 
of reactor and accelerator experiments (see Sec l3.2[) . Solar and atmospheric experiments also exclude 
similar regions, reported in [I]. 

we do not plot it. The determination of primordial Deuterium in eq. ()27j) has to be improved by 
at least a factor of 3 in the uncertainty for it to become competitive in this respect. 

Focusing first on the BBN bounds, the contours for L u = have perfect agreement with 
those in the literature [3J EI] , apart from the fact that we used the most recent cosmological 
values. Introducing the primordial asymmetry L v causes the contours to move upwards and to 
the right, as expected, thus reopening a portion of the parameter space. The bounds in the es 
mixing case for large asymmetries reach lower values of Am 2 4 than those in the fis cases. This 
is because the depletion of v e due to its oscillating into v s after decoupling from the background 
plasma increases the 4 He yield, whereas the corresponding depletion of has no effect on 4 He 
production, since does not interact directly with the protons and neutrons. 

The bounds from LSS are also shifted by the presence of the primordial asymmetry. For a 
given Amy and a given 9 S , the lower final abundances of sterile neutrinos brought about by the 
non-zero asymmetry allows for a larger sterile neutrino mass for a fixed Q u h 2 , resulting in a less 
stringent constraint. 
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Figure 2: The allowed LSND region at 99% C.L. (yellow/light shaded area) compared to the cosmo- 
logical bounds from BBN and LSS in the presence of primordial asymmetries. The darker shaded area 
is already excluded at 99% C.L. by other experiments. The regions below and to the left of the thin 
lines are allowed by BBN because they correspond to Y p < 0.258. The regions below and to the left of 
the thick lines are allowed by LSS because they correspond to Vt v h 2 < 0.8 10~ 2 . 



3.2 Rescuing LSND 

The LSND experiment reported a signal [3] for — > v e oscillations in the appearance of v e in an 
originally beam. The best fit point is located at sin 2 6*lsnd = 3- 1CT 3 , AmJ SND = 1.2 eV 2 (the 
whole allowed region is represented in figJ2]as a lightly colored area). The "3+1 sterile" neutrino 
explanation assumes that the — > v e oscillation proceeds through — > v s — >• u e , therefore 
introducing the large LSND mass scale Aml SND as the mass squared difference of the active 
neutrinos with the mainly sterile state. The effective LSND mixing angle can be expressed in 
terms of the active/sterile mixing angles 9 es and 9^ as 

sin 2 2# LS nd = - sin 2 26 es sin 2 2^ s , (31) 

which reduces to 6*lsnd ~ es ■ 9^ s for small mixing. Both angles are constrained by several other 
experiments that found no evidence of electron or muon neutrino disappearance. As a result [3*%] . 
a large portion of the area indicated by LSND, but not all of it, is ruled out (the darker shaded 
area in fig. |2J). 

The constraint by cosmology are represented by lines in fig. [21 They have been computed 
in the most conservative way, i.e. by choosing for each value of sin 2 2#lsnd the combination 
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1/2 

#es ~ 6 /is ~ ^lsnd ' which is found to minimize the overproduction of sterile neutrinos. When 
these values are ruled out by neutrino data, we set one of the two angles to its largest allowed 
value and determine the other using eq. (J3TJ). On physical grounds, setting 9 es ~ 9^ s means the 
es and fis " channels" both contribute equally to the production of 4 He and Deuterium. From 
eq. (|31j). this implies that the LSND plot in fig. (J2J) can be obtained, in a qualitative sense, from 
fig. © by a re-scaling sin 2 29 — > 1/4 (sin 2 Iff) 2 (for small mixing). This guideline is expected to 
hold for other mixing patterns that involve two flavors with similar mixing angles. 

We plot the isocontours of Y p = 0.258 and Q u h 2 = 0.8 10~ 2 , so that the regions above and 
to the right of these lines are the excluded portions. One recovers the well known result that 
the region excluded by standard BBN (L u = 0, solid thin line) includes all of the allowed LSND 
region. The bound from LSS (solid thick line), updated by the recent cosmological results as 
discussed in Sec J2.31 now also lies below the LSND allowed portions. 

The effect of introducing the primordial asymmetry is represented by the dashed lines: the 
larger the value of the (negative) asymmetry, the more inhibited the production of sterile neu- 
trinos, and the further upper right the 4 He and Q u h 2 contours are displaced. We find that a 
primordial asymmetry L v of the order of — 10 -4 is enough to re-allow the portion of the LSND 
region around (Am 2 SND , sin 2 2#LSND) — (1 eV 2 ,few 10~ 3 ). This result is less stringent than the 
estimate of L Ve ~ O(10~ L - 2 ) in [Hj, that was derived without the use of the explicit neutrino 
evolution equations. 

4 Conclusions 

We investigated the extent to which introducing a primordial leptonic asymmetry {n v ^ n D ), 
much larger than the baryonic one, relaxes the bounds on light sterile neutrinos imposed by 
Big Bang Nucleosynthesis and Large Scale Structure, by suppressing via matter effects the pro- 
duction of the sterile states in the Early Universe. Our momentum-averaged neutrino evolution 
equations took into account the background cosmological evolution and included all the relevant 
features of neutrino mixing, both active-active and active-sterile, in the presence of a primordial 
leptonic asymmetry. 

Working in the framework of constant leptonic asymmetry, we solved numerically the neutrino 
densities as a function of temperature and computed the relevant cosmological observables such 
as the yield of 4 He and Deuterium and the neutrino contribution to the energy density in the 
Universe Q u h 2 , comparing them with the most recent determinations. 

We applied this formalism to the study of the generic cases of a sterile neutrino mixing 
with a pure electron neutrino and with a pure muon one. For different values of the primordial 
asymmetry L u , taken to be the same for all the active flavors, we found the relaxed constraints 
on the parameter spaces. The results are presented in Fig. ^ across a broad range of mixing 
parameters. 

We then studied the specific case of the LSND sterile neutrino, which mixes with a specific 
combination of electron and muon neutrinos. The constraint from standard cosmology in the 
absence of asymmetries are confirmed to rule out the entire LSND region, by both BBN 4 He 
overproduction and the cosmological neutrino mass bound. We found that a primordial asymme- 
try of order C(10~ 4 ) allows the LSND regions around (Am 2 SND , sin 2 2# L snd) — (1 e V 2 , few 10~ 3 ) 
to be reopened. 
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A Appendix 

In this Appendix we lay down explicitly the formalism that we use for the evolution of the 
temperature in our computation. In full generality, a standard analysis [2H EH] gives for the 
differential change of temperature 

dT fdT\ . ( dT\ * 

hH a + hrrr ) N » ( 32 ) 



dt \da J N \dN u 

where the first term accounts for the expansion of the Universe (a being the scale factor) and 
the second term for the variation in the count of the neutrino total degrees of freedom N v . The 
first term can be computed exploiting the conservation of the entropy density in the system of 
photons and e ± per comoving volume 

CC~* 4:7T 2 / Tfl \ 

s =- (p 7 + Pe± + P 7 + P e± ) = (aT) 3 —S [^) , (33a) 

45 /' , J <r \ i 



with S(x) = 1 + — / dy y* UV + ^ 2 + - 7= 5 —7= (33b) 

2vr 4 J y Sy/y^ + x 2 J e Vv 2 +x 2 _|_ i 

that effectively provides ana<->T relation. One obtains (here S = S(m e /T)) 

(9T\ . a ( -TS 1 '' 3 \ Ir , . ( -TS 1 ' 3 \ 

{^) N ^a{j(fsyi-)) = H ^{j(fs^) (34) 

In the same way, one obtains the relation between T v and T 

T,A 3 4 



— = —S. (35) 
T J 11 V ; 

S(m e /T) — 1 effectively parameterizes the entropy in electrons and positrons at any given T, 
including during their annihilation into photons; for T 3> m e , iS — > 11/4. 

The second term in eq. (|32|) can be computed exploiting the conservation of the total energy 
density (at fixed a) 



Ptot = (p 7 + + Pe ± ) = ^T 4 ^2 + (^) 4 + Tl (^) j , (36a) 

60 /* 1 

with 7£(x) = — dy y 2 y / y 2 + x 2 — == (36b) 

7T J e y/y 2 +x 2 _|_ ]_ 

that effectively provides an N v «-> T relation. One obtains the lengthy expression (here 7£ = 
K{m e /T)) 



4/3 



<9T\ • N V T f „ r T dS\ 8 /11\ 4/ * 4/3 ^ r <9ft\ 



(37) 



In our case, this contribution, proportional to JV„, is subdominant with respect to the first one 
in eq. (|32j). since the differential change in N u due to the production of the fourth family of 
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neutrinos is limited by their contribution to the total energy budget of the universe, which is at 
most 0(7/50). This effect is also already included in the expression for the Hubble parameter 
if (jp tot ). Therefore in the n/p evolution equation eq. (J22J) we only keep the first addend of the 
dT/dt relation eq. Finally, we note that in the neutrino evolution equations eq.s (fTlj) we 

have made use of the further approximated expression dT/dt = —H(p tot )T. Eq.s and (J23J) 
reduce to this form at T ^> m e when S reduces to the constant value 11/4. In this same regime 
1Z — > 7/2 and of course T u = T. 
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